clear all;close all; clc; 
R=1.25;  % сопротивление для диода 2Д103А из примера
C=33e-12; % емкость конденсатора задана по условию
tau=R*C;
fprintf('tau=%f ns\n', tau*1e9); % вместо %f записывается то, что стоит после запятой
T=0.01*tau; % Шаг
t=0:T:3*tau; % Создать массив с шагом T от 0 до 3 tau
U(1)=0;
S=10;
alpha=tau/T;
for k=2:length(t);
    U(k)=(S+alpha*U(k-1))/(1+alpha);
end;
plot(t,U);